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Fragmentation functions for hadrons composed of heavy quarks are calculated directly 
rS I from the definitions given by Collins and Soper and are compared with those calculated in 

another way. A new fragmentaion function for a P-wave meson is also obtained and the 
singularity arising at the leading order is discussed. 



Recently it has been shown[l,2,3,4] that the fragmentation functions for a hadron as 
a bound-state of two heavy quarks can be calculated perturbatively. The reason such cal- 
culations are possible is that the bound-states can be described well by the nonrelativstic 
wavefunctions and hence the effect of long distance in the fragmentaion can be factorized 
into the wavefunction, expressed through the radial wavefunction at the origin. Although 
it is realized [2,4] that the fragmentaion functions are universal, the fragmentation func- 
tions calculated in these works are extracted from specific processes. It should be noted 
that in QCD the definitions of the fragmentaion functions (or decay functions) are already 
given by Collins and Soper[5] and these functions should be universal and independent of 
pocesses. The so called factorization thorem is based on such definitions (for the factoriza- 
tion theorem see [6] and references cited there) . It is interesting to calculate fragmentation 
functions directly from these defintions and to compare those obtained previously. On the 
other hand, it is also convenient to work with the defintions for calculating the high order 
corrections. In this letter we will calculate fragmentation functions from the definitions. 

To give the definitions for a fragmentaion function it is convenient to work in the 
light-cone coordinate system. In this coordinate system a 4-vector p is expressed as p^ = 
(p"*",!?", Pt), with p+ = {p^ + p^)/\^, p~ = {p^ —p^)/\/2. Introducing a vector n with 
n^ = (0, 1, Ot), the fragmentation functions for a spinless hadron H are defined as[5]: 

DH/qiz) =|-y"(/a;-e-^^^""/^^Treoior^TrDirac{n-7 < 0|Q(0) 

/•OO 

Pexp{-igs / rfAn-G^(An^)}a;J^(P+,OT)aH(i'+,OT) 
Jo 

/•OO 

Pex.p{igs / dXn ■ G^{Xn^')}Q{0, x", Ot)|0 > 
Dh/g{z) =^^/^^"^"'''^""^' < 0|G''+^(0) 

/•CXD 

{Pexp{-igs / dXn ■ G'(An'^)}}''^al,(P+, OT)aH{P+, Ot) 
Jo 

{Pexp{i^, / rfAn-G'(An^)}}^^G"^'+(0,a;-,OT)|0> 

J x~ 

Where G^{x) = G'^(x)^, G'^(x) is the gluon field and the A"(a = 1, . . . , 8) are the Gell- 
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Mann matrices. The subscription T denotes the transpose. G"''^'^ is the gluon field strength 
and a\j(P) is the creation operator for the hadron H. The function Dh/q{z) or Dh/g{^) 
are interpreted as the probabhty of a quark Q or a gluon G with momentum k to decay 
into the hadron H with momentum component P~^ = zk^ , both are gauge invariant from 
the definitions. The defintions in (1) are the unrenormalized versions, i.e. all quantities 
are bar quantities (for renormalization see [5]). When one has a parton, i.e. quark or 
gluon, instead of the hadron in (1), the perturbative expansion in g^, the strong coupling 
constant, can also be expressed with Feymann diagrams, the Feymann rules can be found 
in [5,6]. The definitions can easily be generalized to hadrons with nonzero spin. 

From the defintions a direct calculation is possible if one can express the hadron 
operators with the parton operators. This is possible for a hadron composed of two heavy 
quarks, where, as mentioned, the nonrelativstic approach works. To see how to express 
the hadron operator as a parton one, let us consider a spinless hadron to be a bound-state 
of a heavy quark Qi and an anti-quark Q2 and to have the quantum number ^Sq- In its 
rest frame the state can be expressed as: 

color 

with 

In Eq(2), al_^{bl^) is the creation operator for Qi{Q2), the indices si and S2 are spin- 
indices, and Aq is a suitable normalization constant. The function ^o(q) = /o(|q|)^oo(^7 '/') 
is obtained through Fourier transformation of the nonrelativstic wavefunction. /o(|q|) goes 
rapidly to zero, when |q| -^ 00. This means, the region with |q| ?^ is dominant in the 
integral of Eq.(2). In this work we will always take nonzero leading order in |q| and the 
hadron mass M = rrii -\- 1712- Through a Lorenz boost the state can be transformed into a 
P ^ state: 

\'So,F>= ^0/ (l^/odqD-^Xs,.,-^ J2 <(Pi)^L(P2)|0 > (4) 

color 
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where 



Pi = : P + Qp, P2 = : P-Qp (5) 

nil + '>T^2 f^i + 1712 



Qp is related through the Lorenz boost to q in the rest frame. Now it is easy to show 
that the bound-state satisfies the normahzation condition given in [5] as required for the 
definitions in (1), if one takes the following form for the hadron creation operator: 
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/rJ n 1 1 

7^/o(|q|)-/|^Xsi.,-^ J2 <(pi)^l(P2) 

color 

1 _ / 2mim2 _ PO + P3 
° "y mi +m2 ' P° + |P| 

Similarly, one can also obtain the operator for the ^Si hadronic state: 

color 

Here a"i(i = 1,2,3) are the Pauli matrices and ei{X){i = 1,2,3) is the component of the 
polarization vector of the hadron in its rest frame. A labels the heilicity with respect to 
the P direction. With Eq.(6) and Eq.(7) we can now calculate the fragmentaion function 
for hadrons with quantum number ^5*0 and ^5*1 from the definions in Eq.(l). In the 
leading order of Qs there are four Feymann diagrams for a quark Qi to decay to those 
hadrons. The diagrams are given in Fig.l, where the double lines represent the Wilson 
line operators in Eq.(l). One can work in the light cone gauge n ■ G = 0. In this gauge 
the Wilson line operators in the defintions Eq. (1) disappear and only one of the four 
diagrams contributes. We will take Feymann gauge. With the Feymann rule in [5,6] one 
obtains the contribution from each diagram to the fragmentaion function Dh/q^{z). It 
is straightforward to calculate, and in the calculation one meets the so called spin-sums, 
which are: 

U{pi, Si)v{p2, S2)Xs,S2 = /^\^^ 75(7 -P -M), 

for H = ^Sq state 

v/^^W^ (8) 

u{pi, si)v{p2, S2){e^{\)a^x)s^s2 = TfTT^'^^ ' ^ ^ ^)^pW ' 7, 

for H = ^Si state 



Here ep(A) is the 4- vector for the polarization of H in the P 7^ frame. Similar results 
for the spin sums in Eq.(8) and below can also be found in [7]. After a straightforward 
calculation we obtain for the H = ^Sq bound-state: 

2 |i?o(0)P 2^2^(1-^)2 



Dh/qA^) 



'81n ml ' (l-yiz)^ 
■ {6 + 18(1 - 2yi)z + {68yf - 62yi + 15)z'^ (9) 



+ 2yi{-18yl + 17yi - 5)z^ + 3yl{l - 2y^ + 2yl)z^} 
and for H = ^Si bound-state: 

Dh/q, [z, a = 0) =^T^ _3 (^s n _....^6 



SItt ml ' [l-yizY 
{2 - 2(1 + 2yi)z + (15 - 22yi + lQyl)z^ 



+ 2yi(-5 + 7yi - Qyl)z^ + 3y?(l - 2y^ + 2yl)z^} 

■ {2 - 2(1 + 2yi)z + (15 - 16yi + 10yl)z^ 

+ 2yi{yi-5)z^ + 3yfz^} 
Here y^ = m,i/M for i = 1, 2 and -Ro(O) is the radial wavefunction at the origin. Comparing 
the fragmentation functions calculated in previous works[3,4] we find agreement. Exchang- 
ing yi -^ y2 one can get D^jiq^. If we take Q2 = Qi, then we obtain the fragmentation 
functions for the quarkonia (QiQi). 

Now we turn to gluon fragmentaion function. If a hadron has Q2 = Qi, then a gluon 
can decay into that hadron at the order of g^. In this leading order there are also four 
Feymann diagrams depicted in Fig. 2. Since Hadrons are colorless, there are always two 
gluons attached to the fermion lines to bulid a color singlet. Because charge- conjungat ion 
invariance a gluon can not decay at this order into a hadron with the quantum number 
"^5'i(i.e. J^*^ = 1~~). Computing the contributions from the diagrams we obtain the 
fragmentaion function for a ^5*0 hadron: 

Dh/g{z) = ^^^^{(3 - ^z)z + 2(1 - ^)ln(l - z)} (11) 



This is also in agreement with that obtained in [2]. With our results calculated directly 
from Eq.(l) one can start with the defintions to calculate the high order corrections. It 
should be pointed out that the fragmentaion functions calculated here are at the scale 
H = M. For the functions at an arbitrary scale /i one needs to solve the correponding 
renormalization group equations and use the results in Eq.(9,10,ll) as initial conditions. 

With Eq.(9,10,ll) we complete the calculations for the fragmentation functions of a 
5"- wave meson. The long-distance effect in these functions are included in the nonrelativstic 
approach through the radial wave functions at the origin, which can be calculated from 
potential models or can be directly measured in leptonic decays. For a P-wave meson the 
calculation is lengthy but still straightforward. Considering a meson to be a J^'^ = l"*""*" 
bound state of (QiQi), the creation operator for this hadron is: 



a^(P,A) = Ao^|^y"^/i(|q|)^e.,(A)Kx).,..£n...3<(Pi)^>lfe) (12) 



and the spin sum used in the calculation: 



w(Pi, si)v{p2, S2)qi^ei^{X){cTi.,x)s^s2£ii 



*2«3 



^ 4^(7 ■ Pi + m^)e^,^pP'^e;{X)q;Yil ■ P - M){j ■ p^ - rm) 

At order a^ a gluon can decay into the J^*^ = 1 "'""'" meson, the Feymann diagrams are 
same as in Fig. 2. Keeping the leading order at |q| we obtain: 

^ (14) 

DnM^, A = ±1) = Y^^{ (Y37) (1 - + -^)} 

The -Ri(O) is the derivative of the radial wavefunction at the origin. The results in Eq.(14) 
are divergent when z -^ 1. The reason is that the gluon exchanged between the quarks 
in Fig. 2 can move collinearly along the quarks, if we take the nonzero leading order at 
|q|. Such singularity is actually well known as the singularity in the zero-binding limit [8], 
appearing in the calculations of the hadronic decay widths of a P-wave meson. The situ- 
ation is similar, if we transpose the diagrams in Fig. 2 and calculate the gluon distribution 



of the hadron. The singularity at z = 1 means that the long-distance effect is present not 
only in the -Ri(O) but also in the remaining parts which we calculated perturbativly and 
hence it prevents the naive use of perturbation thoery for z near to 1. Recent progress [9] 
on the P-wave hadron decay shows that the long-distance effect can be clearly seperated 
from the short-distance effect by a new factorization theorem[9] for the decay, where one 
realizes that any meson is a superpostion of many components: 

\M >= vCqqIQQ > +^QQa\QQG >+■■■ (15) 

and the configuration of \QQG > plays an important roll in the decay. In the hadronic 
decay of a P-wave meson one should also take account on the \QQG > configuration and 
the decay width can be written in the new factorization form, where instead of one, two 
parameters represent the long-distance effect and the remaining parts can be calculated 
perturbativly without the singularity. To apply this idea in our case, one needs to know 
the wavefunction "^qQq in Eq.(15), which can not be obtained through the nonrelativstic 
approach. A detailed study is needed. Here we recall a phenomenological treatment for the 
singularity at z = 1. Similar treatment of the decay was also used in [8,10]. In calculating 
the contribution from the diagrams to Du/g w6 integrated the transverse monmentum 
|Pt|, which is carried by the exchanged gluon, from zero to infinity. However, the meson 
is not a point-like particle and it has a spatial extension. The extension in the transverse 
direction is of the order R ^ 1/M. If |px| is smaller than 1/-R, one should treat the gluon 
as a part of the meson. From this point of view, the integral over |pt| should be from a 
nonzero Ipxlmin ~ 1/-R to oo. We introduce a parameter 13: 

|PT|min=/3M (16) 

If the picture given above is correct, 13 should be around 1. Here we take it as a free 



parameter. With Eq. (16) we obtain: 



I'«/ofeA = 0) = i|^iM!((l_,)2 + ,2^.)-3 



108 4nmi 

■ {z{-2z'^ + \\z^ - 26z^ + 35z^ - 30z^ + ITz^ -Qz + l) 

+ W'^z^{-2z^ + 7z^ - 9^2 + 5z - 1) + ?>[5^z^{-2z^ + Az^ - 3z + 1)} 

■ {z{-z'^ + Qz^ - IQz^ + 252^ - 2hz^ + IQz^ -6^ + 1) 

+ 3/3^^3 (-^^ + Az^ - 7z^ + Iz^ -4^ + 1)} 

(17) 

The results in Eq. (17) are finite at ^ = 1 and for z ^ the expressions approach those 

given in Eq. (14). In Eq.(17) there are two parameters which need to be determined, the 

i?']^(0) can be obtained from potential models or measured elsewhere, while /3 is unknown 

in principle. However, in the fragmentation functions the dependence on the momentum 

fraction z is predicted by the perturbative calculations. It should be stressed that it 

is important to know the fragmentation functions for the P-wave meson or quarkonia. 

If we take Q\ as a c quark, the quarkonia then corresponds to a Xci meson and Xci 

can substantially decay into J j'^ via radiative decay. Hence, a theoretical prediction 

of Xci production via fragmentation is very important for the prediction for the J j'^ 

production(for example see [11] and the references cited there). 

The defintions given in Eq.(l) are for unpolarized partons. One can generalize the 
definitions to polarized partons and then obtain fragmentation functions for a polarized 
parton. In this way one can study through heavy meson production the polarization 
properties of the quarks and gluons produced in high energy collisions, where one expects 
for large Pt processes that the heavy mesons are produced dominantly via fragmentation[2]. 

To summarize: We have shown in this letter that the fragmentation functions for a 5"- 
wave meson composed of two heavy quarks can be calculated directly from the definitions 
given by Collins and Soper in [5] and the results obtained are same as those obtained by 
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other means. A new fragmentation function for a P-wave meson is also obtained, but in 
contrast to a 5'-wave meson, there is an extra parameter which needs to be estimated. 
From the naive picture given above this parameter should be around 1. Finally within 
the defintions one can conveniently study higher order corrections to the fragmentaion 
functions. 
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Figure Captions 

Fig.l. The Feymann diagrams for a quark to decay into a hadron. 
Fig. 2. The Feymann diagrams for a gluon to decay into a hadron. 
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Fig. 1 
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Fig. 2 
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